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Abstract 



The prospects for determining \Vub\ from exclusive B semileptonic de- 
cay are discussed. The double ratio of form factors (/(^~*'')//(-^^^*))/ 
^j(D^p) I f{D^K')^^ is calculated using chiral perturbation theory. Its devi- 
ation from unity due to contributions that are non-analytic in the symmetry 
breaking parameters is very small. Combining experimental data obtainable 
from B pii'i, B — > K*£i and D pi ft, can lead to a model independent 
determination of \ Vub\ with an uncertainty from theory of about 10%. 
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The next generation of B decay experiments will test the flavor sector of the standard 
model at high precision. The basic approach is to determine the elements of the CKM 
matrix using different methods and then check for the consistency of these results. At 
the present time CP non-conservation has only been observed in kaon decay arising from 
j^o _ J^o ]2iixing. Many extensions of the minimal standard model (e.g., models with several 
Higgs doublets or low energy supersymmetry) have new particles with weak scale masses 
that contribute to flavor changing neutral current processes like K'^ — mixing, 
mixing, B K*'y, etc., at a level comparable to the standard model. 

At the present time, the magnitude of the h u CKM matrix element is determined by 
comparing experimental results on the inclusive electron spectrum in the endpoint region 
with phenomenological models [||], or by comparing experimental results on B ^ p i Uf and 
B ^ IT £ i?i with phenomenological models and lattice QCD results 0. These two approaches 
yield remarkably consistent determinations of \Vub\, but have large uncertainties. 

This paper is a comment on the proposal to determine \Vub\ using a combination 
of heavy quark symmetry and SU{3) flavor symmetry. The basic idea is to compare 
D K* iui with the Cabibbo suppressed decay D — > plvf. Using heavy quark symmetry 
the S'f/(3) violations in the form factors that occur in these decays are related to those 
that occur in a comparison of -B ^ K*ii (or B K* UfUf) with B pii^e- Therefore, 
experimental data on B ^ K*ii in conjunction with data on D ^ piui and D -—>■ K* lu^ 
can be used to determine \Vub\- This proposal is complementary to other approaches for 
determining \Vuh\i since it relies on the standard model correctly describing the rare flavor 
changing neutral current process B K*ii. 

In this letter we compute corrections to these form factor relations which violate both 
chiral and heavy quark symmetry, and are non-analytic in the symmetry breaking param- 
eters. We also reconsider the influence of long distance effects on the extraction of the 
B K* form factors from B —>■ K*i£. 

We denote the form factors relevant for semileptonic transitions between a pseudoscalar 
meson P^'^\ containing a heavy quark Q, and a member of the lowest lying multiplet of 
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vector mesons, V, by g^^^^\ j{H-*v) a^^^^\ where 

{V{p\ e) I Q \H{p)) = i ^^^^^^^ £,.a. e*^ (p + p')' {p - p'Y , (1) 
(y{p\ e) I g7,75 Q \H{p)) = /(^^^^ e; + (e* ■ p) (p + p'), + aL^"^^^ (e* ■ p) {p - p'), , 

and _ —^^-^23 = 1. We view the form factors f and a± as functions of the dimension- 
less variable y = v-v', where p = tuh v, p' = my v', and = {p—p'Y = 'm'jj+my—2'mH rny y- 
(Although we are using the variable v ■ v', we are not treating the quarks in V as heavy.) 
The experimental values for the D K* i z/£ form factors assuming nearest pole dominance 
for the dependences are 

'1.9±0.1)GeV 



1 + 0.63 (y-l) 



{D^K^f ^ (0.18 ± 0.03) GeV'^ 
~ 1 + 0.63 {y-1) ' 

The shapes of these form factors are beginning to be probed experimentally The form 
factor a_ is not measured because its contribution to the D —>■ K* i ui decay amplitude is 
suppressed by the lepton mass. The minimal value of y is unity (corresponding to the zero 
recoil point) and the maximum value of y is (m|, + m\,) / {2m£, Tn^*) — 1-3 (corresponding 
to = 0). Note that f{y) changes by less than 20% over the whole kinematic range 
1 < y < 1.3. In the following analysis we will extrapolate the measured form factors to the 
larger region 1 < y < 1.5. The full kinematic region for B— * pi ij£ is l<y<3. 5. 

The differential decay rate for semileptonic B decay (neglecting the lepton mass, and 
not summing over the lepton type i) is 

dT{B^piiy,) Gl\Vut\' 



dy 48 7r3 

Here S^^^^\y) is the function 



■niB 



mlS^^-^\y). (3) 
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f 



(y) + - 6yr + 3r^ 



+ 4Re[af r {y - r){y' - 1) 



+ 4 



(y) 'm%r^{y'^ 



(y) m%r\l + r'-2yr){y'-l) 



(2 + y2 - 6yr + 3r') [1 + 



with r = mv / rnu- The function 5^^^^^ depends on the ratios of form factors 



(4) 



and g'y^^^) I f^^^^\ S'^^^'^\y) can be estimated using combinations of SU{'i) flavor sym- 
metry and heavy quark symmetry. SU{'i) symmetry imphes that the p+ form factors 
are equal to the B K* form factors and the B^ 

times the B — > K* form factors. Heavy quark symmetry implies the relations [| 



p° form factors are equal to l/\/2 
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(y) 



TTLd) 



rriB, 



niB 



as 


[rrib) 




[rric) 


as 


{nib) 


as 


[rric) 


as 


[rrib) 


as 


[mc) 



-6/25 



/ 



-1 -6/25 



a 



(D^K*) 



-1 -6/25 



The second relation is obtained using a'^^^ ■* 
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(y) 



(5) 



— a+ ', valid in the large rric limit. 
Using Eq. (|^) and SU{3) to get B^ ii^e form factors (in the region 1 < y < 1.5) 

from those for D K*lv£ given in Eq. (0) yields S'^^^P^y) plotted in Fig. 1 of Ref. §. 
The numerical values in Eq. @) differ slightly from those used in Ref. |^ . This makes only a 
small difference in S^^^p\ but changes 5^^"*'''' more significantly. In the region 1 < ?/ < 1.5, 



\5'^^^P\y)\ defined in Eq. (^) is less than 0.06, indicating that 



a 



and g^^^P^ make a 



small contribution to the differential rate in this region. 

This prediction for S^^'^p^ can be used to determine \Vub\ from a measurement of the 
B pd-Vi semileptonic decay rate in the region 1 < y < 1.5. This method is model 
independent, but cannot be expected to yield a very accurate value of \Vub\- Typical SU{?>) 
violations are at the 10 — 20% level and one expects similar violations of heavy quark 
symmetry. 

Ref. proposed a method for getting a value of S^^~*p^ (y) with small theoretical uncer- 
tainty. They noted that the "Grinstein-type" double ratio 



R{y) = (6) 

is unity in the limit of SU{3) symmetry or in the hmit of heavy quark symmetry. Corrections 
to the prediction R{y) = 1 are suppressed by ms/nic^h {jnu,d ^ "^s) instead of ms/AgcD or 
Aqcd/'"^c,6- Since R{y) is very close to unity, the relation 



mB - rriK* 



(7) 



together with measurements of |/^^^'^*'*|, \f^^^''^\, and S^^^^*^ will determine S^^^^^ with 
small theoretical uncertainty. The last term on the right-hand-side makes Eq. equivalent 
to Eq. (P) in the y — > 1 limit. The ratio of the (2 + — 6yr + 3r^) [1 + S^^^^\y)] terms 
makes only a small and almost y-independent contribution to S^^^^^ / S^^^-^*^ in the range 
1 < y < 1.5. Therefore, corrections to Eq. are at most a few percent larger than those 
to R{y) = 1. 

j^g^g oiYQg^dy been determined. \f^^^^^\ may be obtainable in the future, for 
example from experiments at B factories, where improvements in particle identification 
help reduce the background from the Cabibbo allowed decay. The measurement B{D — >■ 
pHvi)/B{D K*Hue) = 0.047 ±0.013 |[ already suggests that |/{^-p)//(^-^*)| is dose 
to unity. Assuming SU{3) symmetry for the form factors, but keeping the explicit my- 
dependence in S^^^^\y) and in the limits of the y integration, the measured form factors 
in Eq. (|) imply B{D p^Ivt)/B{D K*^Ivi) = 0.044.[| is obtainable from 

experimental data on B K* oi B ^ K*ii. While the former process is very clean 
theoretically, it is very difficult experimentally. A more realistic goal is to use B — > K*i£, 
since CDF expects to observe 400 — 1100 events in the Tevatron run II (if the branching 
ratio is in the standard model range) 0. There are some uncertainties associated with long 
distance nonperturbative strong interaction physics in this extraction of S^^^^*\y). To use 



^This prediction would be |K:d/V^sP/2 ~ 0.026 with nip = nix*- Phase space enhances D ^ p 
compared to D ^ K* to yield the quoted prediction. 
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the kinematic region 1 < y < 1.5, the form factor ratio f^^^f^f^^^^*) in Eq. must be 
extrapolated to a greater region than what can be probed experimentally. For this ratio, 
the uncertainty related to this extrapolation is likely to be small. 

The main purpose of this comment is to examine the deviation of R from unity using 
chiral perturbation theory. We find that it is at the few percent level. The uncertainty 
from long distance physics in the extraction of S^^^^*^ is also reviewed. On average, in the 
region 1 < y < 1.5, this is probably less than a 10% effect on the B K*ii decay rate. 
Consequently, a determination of \Vub\ from experimental data on D ^ K*ii'i, D pii^e, 
B K*ii and B ^ pii?e with an uncertainty from theory of about 10% is feasible. 

(i) Chiral perturbation theory for R 

The leading deviation of R from unity can be calculated using a combination of heavy 
hadron chiral perturbation theory for the mesons containing a heavy quark and for the lowest 
lying vector mesons. We adopt the notations and conventions of Refs. |TO|JTT[] . The weak 



current transforms as (3l, and at the zero recoil kinematic point there are two operators 
that are relevant for P^'^) y transition matrix elements (where P*-^^ = B, P^'^^ = D, and V 
is one of the lowest lying vector mesons p, u, K*, 0). Demanding that the Zweig suppressed 
Dg —>■ uj iui process vanishes relates the two operators, yielding |T^ 

Qa 7m(1 - 75) Q = Tr[^],7M(l - 75)^^^'^^^] • (8) 
Here repeated SU{3) indices are summed and the trace is over Lorentz indices. contains 



the ground state heavy meson doublet, is the nonet vector meson matrix ||rT[, and /5 is a 
constant. The leading contribution to -R(l) — 1 arises from the Feynman diagrams in Fig. 1. 
Diagrams with a virtual kaon cancel in the double ratio R. Neglecting the vector meson 
widths,^ these diagrams yield 



^The only significant width is that of the p meson. Since it occurs in the loop graph involving an 
y, neglecting the p width amounts to treating T p/2mr^ <^ 1, which is a reasonable approximation. 
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FIG. 1. Feynman diagram that gives the leading contribution to R{1) — 1. The dashed hne is 
a TT or an r]. The black square indicates insertion of the weak current. 



i?(l)-l = - 



12 7r2 /2 

where A^^^^ = uib* — itib, A^^^ = m^* — uid, and for m > A, 



G{m,A) = Y^-- -^—^ — arctan 1^ — j-A^lnm. (10) 

Here g2 is the ptuvr coupling, g is the DD*tt coupling, and / ~ 131 MeV is the pion decay 
constant. In the nonrelativistic constituent quark model g = g2 = ^ while in the chiral 
quark model [0 g = g2 = 0.75. Experimental data on r ^ u; tt z/,- in the region of low un 
invariant mass gives g2 — 0.6 ||14| . 

For small A, Eq. for R{1) — 1 has a non- analytic y/mg dependence on the light 
quark masses. This cannot arise from corrections to the current in Eq. (H) or to the chiral 
Lagrangian, and must come from 1-loop diagrams involving the pseudo-Goldstone bosons 
TT, K, 1]. Using the measured values of the pion and eta masses, Eqs. (H) and ([T0|) imply 
R{1) = 1 — 0.035(7(72- There may be significant corrections from higher orders in chiral 
perturbation theory. However, the smallness of our result lends support to the expectation 
that -R(l) — 1 is very close to zero. There is no reason to expect any different conclusion 
over the kinematic range 1 < y < 1.5. 

(ii) Long distance effects and the extraction of from B — > K*ii 

The decay rate for B K* vn vn could determine S^^^^*^ free of theoretical uncer- 
tainties. However, experimental study of this decay is very challenging. A more practical 
approach to extracting this quantity is to use B K*ii. The differential decay rate is 



: — - 1 I tils '/t^^* 



dy 24 7r3 V47ry 



1^9(2/) P+|CioP [1 + A 



5(^-^^*)(y)[l + %)], (11) 



This and Eq. (0) allow us to rewrite Eq. as 



dl/ \V,lVt,\'' «2 |^g^^)|2 + |c^^|2 1 + A(y) 1 + d{y) 



^ ml / TUB - nip ^"^ 



dy ' ^''^ 



which can be directly used to extract \Vub\- Unitarity of the CKM matrix implies that 
IV^*!^^! — iKs^cbl with less than a 3% uncertainty. The fine structure constant, a = 1/129, 
is evaluated at the VT-boson mass. d{y) parameterizes long distance effects, and will be 
discussed below. A{y) takes into account the contribution of the magnetic moment operator, 
O7 = (e/167r^) nib (sl a^y bn) F^" (a factor of —AGFV^gVtb/ has been extracted out in the 
definition of operator coefficients). Ref. 0] (see also Ref. ||15|) found using heavy quark 
symmetry that A(?/) ~ —0.14 — 0.08(?/ — 1) in the region 1 < y < 1.5. Corrections to 
this are expected to be small since there are no l/m^ corrections to A(l). Cio is the 
Wilson coefficient of the operator Oio = (e^/167r^) (sl 7^ 6^) (17^^75 £). C^iy) takes into 
account the contribution of the four-quark operators, Oi — Og, and the operator O9 = 
(e^/167r^) (sl7^6/,) {i'j^i). In perturbation theory using the next-to-leading logarithmic 
approximation 

C,iy) = C, + hiz, y) (3Ci + C2 + 3C3 + C4 + 3^5 + C,) -\ /i(0, y) (C3 + 3C4) 

- i /i(l, y) (4^3 + 4^4 + 3C5 + Ce) + ^ (3C3 + C^ + 3C, + Cg) , (13) 

where z = rric/nib. Here 



,/ ^ 8, 8 4 2 , , r I In ^ + Vl ^^ _ 

h{u,y) = ~-\nu + — + -x--i2 + x)^J\l-x\l 1 - (14) 

[ 2 arctan(l/Va; — 1) ; a; > 1 , 
where x = Au'^ml/{ni^Q + mj^, — 2mBniK* y)- Using mt = 175 GeV, nib = 4.8 GeV, nic = 
1.4 GeV, as{niw) = 0.12, and as{nib) = 0.22, the numerical values of the Wilson coefficients 
are Ci = -0.26, C2 = 1.11, = 0.01, C4 = -0.03, C5 = 0.008, Cq = -0.03, C7 = -0.32, 
Cg = 4.26, and Cio = —4.62. Of these, Cg and Cio are sensitive to nit (quadratically for 
mt > mw)- 
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In Eq. (0) the second term on the right-hand-side, proportional to h{z, y) comes from 
charm quark loops. Since the kinematic region we are interested in is close to = 4m^, 
a perturbative calculation of the cc loop cannot be trusted. Threshold effects which spoil 
local duality are important. It is these long distance effects that give rise to the major 
theoretical uncertainty in the extraction of \Vub\ from the B — » K*ii differential decay rate 
using Eq. The influence of this long distance physics on the differential decay rate is 

parameterized by d{y) in Eq. ([TTD , where setting d{y) = gives the perturbative result. 

For the part of the c c loop where the charm quarks are not far off-shell, a model for h{z, y) 
which sums over 1 cc resonances is more appropriate than the perturbative calculation. 
Consequently, we model the part of h{z, y) with explicit g^-dependence in Eq. (|1^) with a 
sum over resonances 



17 1 calculated using factorization 



_8 8 37r^ ^ r^wg(^(")^£l) 

The resonances ip^""^ have masses 3.097 GeV, 3.686 GeV, 3.770 GeV, 4.040 GeV, 4.160 GeV, 
and 4.415 GeV, respectively, and their widths r^(n) and leptonic branching ratios B{ip^"^ — >■ 
ii) are known [|1^]. The factor k, = 2.3 takes into account the deviation of the factorization 
model [|19[ parameter 02 from its perturbative value. Denoting the value of Cg{y) in this 
model by Cg{y), its influence on the differential decay rate is given by d{y) defined as 

\C',{y)\' + ICioP = i\C,iy)\' + ICioH [1 + d{y)] . (16) 

d{y) is plotted in Fig. 2 (solid curve). The physical interpretation of the 1 resonances 
above 4 GeV is not completely clear. It might be more appropriate to treat them as DD 
resonances than as cc states. It is possible that for these resonances factorization as modeled 
by Eq. (p!5| ) with n = 2.3 is not a good approximation. Including only the first three 1 
resonances in Eq. (|1^), yields d{y) plotted with the dashed curve in Fig. 2. 



^The four-quark operators involving light n, d, and s quarks also have uncertainty from lon^ 
distance physics. However, this is expected to have a very small effect on the B — > K*ii rate. 
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FIG. 2. d{y) defined in Eq. (|T^). The solid curve takes into account all six 1 cc resonances 
according to Eq. ([l5|), whereas the dashed curve is obtained including only the three lightest ones. 

This estimate of d{y) based on factorization and resonance saturation differs from that 
in Ref. in two respects. Firstly, the phase of k is viewed as fixed because recent data 
has determined the sign of the ratio of factorization model parameters, a2/ai, and the 
phase of ai is expected to be near its perturbative value ||2^. Secondly, since the resonance 
saturation model only represents the c c loop for charm quarks that are not far off-shell, we 
have only used it for the part of h{z,y) in Eq. (|T3p with explicit dependence, retaining 
the perturbative expression for the first two terms, —(8/9) In 2; + 8/27. The Inz term has 
dependence on rrib, which is clearly short distance in origin. This reduces somewhat the 
magnitude of d{y) and makes it more symmetric about zero (compare Fig. 2 with Fig. 6 of 
Ref. 0). It would be interesting to have a more physical separation between the long and 
short distance parts of the amplitude. 

Whether it is reasonable to use factorization for the resonances above 4 GeV can be tested 
experimentally, since these states cause a very distinctive pattern in dT/dy. In Fig. 3 the 
shape of dT/dy is plotted in the region 1 < y < 1.5 using the resonance saturation model 
for d{y) (solid curve). Experimental support for this shape would provide evidence that 
this model correctly describes the long distance physics parameterized by d{y). Although 
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1.1 1.2 1.3 1.4 1.5 

y 



FIG. 3. dT{B K*£i)/dy in units of [\C9{1)\^ + |CioP] m^m^, G||y/,yi6|2aV(3847r5) as 
given in Eq. (pT]). The solid curve takes into account all six 1 cc resonances, the dashed curve 
includes only the three lightest ones, and the dotted curve is the perturbative result (i.e., d{y) = 0). 

d{y) gets as large as ±0.2, since it oscillates, its influence on the B K*i£ decay rate 
in the region 1 < ?/ < 1.5 is about —8% compared to the perturbative result (which is 
plotted with the dotted curve in Fig. 3). Even if our estimates of this long distance physics 
based on factorization and resonance saturation has a 100% uncertainty (a prospect that 
we do not consider particularly unlikely), it will only cause about a 4% uncertainty in this 
determination of \Vub\- Including only the first three 1 resonances in the sum in Eq. (p!5D 
yields the dashed curve in Fig. 3. In this case d{y) causes a —13% change in the B K*ii 
decay rate in the region 1 < y < 1.5. 

(iii) Nearer term prospects 

Without information on the y spectrum for the B decay rates in Eq. (|T^, it is still 
possible to determine \Vub\ by comparing the branching ratios ioi B piui and B K*£i 
in the region 1 < y < 1.5. Integrating Eq. (|T2p over 1 < y < 1.5 we can write 
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y<'-^ \Vt*sVtt\^ «2 cl + c!o (i + A) (i + rf) 



m 

X 



/(15-P)(1) 



T{B ^ K*i 



y<1.5 



Here the barred quantities, Cg, A, and d denote the averages of ICgd/)^, A{y), and 
weighted with S^^^^*\y). Using the shape for ^C-^^-f^*) predicted from heavy quark sym- 
metry, we find Cg = 4.58, A = —0.16, and d = —0.08. Note that the y-dependence of Cg is 
small and Cg is close to Cg. In Eq. ( |17|) the y-dependence of the ratio f'^^^^\y)/f^^^^*\y) 
has been neglected. If the shape of these form factors can be approximated with a pole 
form, then the pole masses of 2.56 GeV for f(D-^K*) 2Ah GeV for f^^~^p^ (correspond- 
ing to D** and to D**) imply that varies by less than 1% over 
the range 1 < y < 1.5. SU{'i) symmetry and the measured D K* form factors imply 
that contributes only about 23% of the D ptvt decay rate. Using this predic- 
tion for and assuming that f^^^p^ and fi^-^^*) have the same ^/-dependence, yields 
B{D pHue)/B{D ^ K*'>£ue) = 0.044 |/(^-^'')(l)//(^^^*)(l)|2. 

In the region = {pi + piY < rrij^^ (corresponding roughly to y > 2), one cannot 
use the double ratio and Eq. (|12[). Moreover, the O7 contribution to the B K*ii rate 
is large and proportional to 1/g^, so the (leading order) heavy quark symmetry relations 
between the tensor and (axial-)vector form factors^ introduce a significant uncertainty. For 

< rrij^^, one can do better using SU{3) flavor symmetry alone to predict dr(_B — >■ 
Tilvij/dq^ from a measurement of dr(i? Kii^jdc^. Since this region is far from g^ax) 



the B* pole contribution ||2^ is unlikely to upset the SU (3) relations. The O7 contribution 
to dr(i? Kl€)/(\q^ is at the 10 — 15% level, fairly independent of q^. In the region 
(1 - 2)GeV^ < < m^/^, neglecting rn\.^/m\, 



dr(50 ^ 7r+£z/£) \VubV Stt^ 1 dV{B^K> 



dg2 \V^%Vtb\'' |C9(g2)+2C7|2 + |Cio|2 dg2 

A similar relation also holds for integrated rates. 



(18) 



It was argued in Ref. |21| that heavy quark symmetry can be used even at small g^. 
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A measurement of the B K*ii decay rate is unlikely before the Tevatron run II. 
Without this measurement, one has to rely on predicting the B —>■ p form factors from 
D —>■ p using heavy quark symmetry, or from D K* using both chiral and heavy quark 
symmetries. As discussed following Eq. (|^, recent experimental data suggests that the 
SU (3) relation between f(D~^K*) y(D-»p) violated by more than 15%. Heavy quark 

symmetry and the measured D K* form factors in Eq. (H) imply that the B^ p^ii^e 
branching ratio in the region 1 < y < 1.5 is 5.9 |Kt6p. The measured decay rate B{B^ — >■ 
p+iui) = (2.5 ± 0.4l[J:^ ± 0.5) X 10~^ together with iKfcl ~ 0.003 imply that about 20% 
of B^ — > p^iue decays are in the range 1 < y < 1.5. 

Despite the presence of long distance effects associated with the cc resonance region, the 
B — > K*£i rate can be used in Eq. ([T^ ) to determine \ Vub\ with a theoretical uncertainty that 
is about 10%. Experimental verification of the distinctive y-dependence of the differential 
rate associated with the 1 resonances above 4 GeV (see Fig. 3) would reduce the theoretical 
uncertainty from long distance effects. A precise value of \Vub\ may be available from other 
processes, e.g., the hadronic invariant mass spectrum in inclusive B — > X„£z/^ decay |^ or 
from lattice QCD results on exclusive form factors before the B — > K*ii decay rate 
is measured. In that case, Eq. (|T2p gives an accurate standard model prediction for the 
B K*ii decay rate in the region 1 < y < 1.5. Comparison with data may signal new 
physics or provide stringent constraints on extensions of the standard model. 
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